function drawDeltas(hP, deltaValues, color, lineW, units)
    if ~isempty(deltaValues)
        xLim = xlim(hP);
        yLim = ylim(hP);
        scale = get(hP, 'XScale');
        width = xLim(2) - xLim(1);
        height = yLim(2) - yLim(1);
        if strcmp(units, 'rad/s')
            deltaValues(:,1) = deltaValues(:,1)*2*pi;
        end

        for i = 1:length(deltaValues(:,1))
            x = [deltaValues(i,1), deltaValues(i,1)];
            y = [deltaValues(i,2)*1000, deltaValues(i,2)*1000 + height * 0.11];
            line(x, y, 'Parent', hP, 'Color', color,  'LineWidth', lineW);
            if strcmp(scale, 'log')
                x = [deltaValues(i,1), deltaValues(i,1)*(10^(0.005*log10(width)))];
            else
                x = [deltaValues(i,1), deltaValues(i,1) - 0.01*width];
            end
            y = [deltaValues(i,2)*1000 + height * 0.11, deltaValues(i,2)*1000 + height * 0.11 - height * 0.033];
            line(x, y, 'Parent', hP, 'Color', color,  'LineWidth', lineW);
            if strcmp(scale, 'log')
                x = [deltaValues(i,1), deltaValues(i,1)/(10^(0.005*log10(width)))];
            else
                x = [deltaValues(i,1), deltaValues(i,1) + 0.01*width];
            end
            line(x, y, 'Parent', hP, 'Color', color,  'LineWidth', lineW);
        end
        xlim(hP, xLim);
        ylim(hP, yLim);
    end
end